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By  deriving  a  dynamical  differential  equation  for  the  electron  distribution  function  in  the  presence  of  a 
nonadiabatic  sequential-tunneling  current  under  an  ac  electric  field  through  a  multiple-quantum-well  system, 
the  self-consistent  Hartree  model  is  generalized  for  the  calculation  of  electronic  states  with  the  inclusion  of 
nonadiabatic  effects  (dependence  on  the  time  derivative  of  the  applied  ac  electric  field)  in  quantum  wells.  The 
influences  of  different  doping  profiles,  temperatures,  and  amplitudes  of  an  applied  ac  electric  field  bn  the 
electron  distribution  function  and  sequential  tunneling  are  studied.  This  work  provides  a  fully  quantum- 
mechanical  explanation  to  the  previously  proposed  current-surge  model  to  a  leading-order  approximation. 
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L  INTRODUCTION 

Recently,  transient  transport  properties  of  semiconductor 
quantum  wells,1'3  superlattices4’5  and  quantum  dots6  have 
been  a  subject  of  interest  Various  nonlinear  properties  hav¬ 
ing  their  origin  in  the  Coulomb  interaction  have  attracted  a 
great  deal  of  attention.3’7'10  Self-sustained  current  oscilla¬ 
tions  and  multistability  have  been  predicted  in  tunneling  cur¬ 
rents  through  doped  semiconductor  superlattices  and  mul¬ 
tiple  quantum  wells  (MQW’s).7  They  are  attributed  to  the 
dynamics  of  domain  walls  separating  the  electric-field  do¬ 
mains.  In  addition,  oscillations  in  the  sequential-tunneling 
current  have  been  predicted  in  MQW’s,  even  in  the  absence 
of  electric-field  domains  due  to  nonadiabatic  effects.10  The 
nonadiabatic  effects  discussed  in  this  paper  are  associated 
with  the  fact  that  a  transient  conduction  current  depends  not 
only  on  an  electric  field  but  also  on  its  time  derivative  due  to 
quantum-well  capacitive  coupling.  The  quantum-well  capaci¬ 
tance  is  of  the  order  of  picofarads.  However,  the  resistance  of 
the  MQW  sample  with  a  thick  barrier  between  the  wells  used 
in  this  study  is  of  the  order  of  teraohms  due  to  the  extremely 
small  sequential-tunneling  current.  As  a  result,  the  nonadia¬ 
batic  effect  occurs  on  a  time  scale  of  seconds,  which  makes 
electron  tunneling  depend  on  the  time  derivative  of  the  ap¬ 
plied  electric  field  in  addition  to  the  field  itself  for  low  ac 
frequencies  of  the  order  of  a  few  hertz.11,12 

In  this  paper,  we  consider  the  sequential- tunneling  trans¬ 
port  of  electrons  in  an  MQW  system  in  the  presence  of  an 
applied  ac  electric  field.  We  assume  that  the  lattice  tempera¬ 
ture  is  kept  constant  and  the  electrons  are  in  thermal  equilib¬ 
rium  with  the  lattice,  so  that  the  electron  temperature  is  the 
same  as  that  of  the  lattice.  The  sequential  tunneling  of  elec¬ 
trons  (of  the  order  of  nanoseconds)  through  a  thick  barrier 
between  adjacent  quantum  wells  in  an  MQW  system  is  a 
very  slow  process  as  compared  to  the  coherent  tunneling  of 
electrons  (of  the  order  of  sub-picoseconds)  through  a  thin 
barrier  in  a  superlattice  system.  However,  electrons  during 
the  sequential-tunneling  process  still  “see”  an  instantaneous 
electric  field  because  of  Tt<27rlFl  for  low  ac  frequency  fl 
(of  the  order  of  a  few  hertz)  with  rt  being  the  sequential- 
tunneling  time  (of  the  order  of  nanoseconds). 
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Adiabatic  electrons  in  an  MQW  system  with  an  applied  ac 
electric  field  stay  in  the  equilibrium  states  with  a  constant 
Fermi  level,  since  re<rti  with  re  being  the  energy- 
relaxation  time  due  to  the  very-long  sequential-tunneling 
time  within  which  an  equilibrium  state  can  be  established  by 
the  much  faster  inelastic  scattering  of  electrons  inside  the 
quantum  well.  However,  the  electron  density  can  vary  with 
time  if  the  electrons  in  the  quantum  well  stay  in  the  nona¬ 
diabatic  state. 10,13,14  For  the  nonadiabatic  state,  the  nonadia¬ 
batic  effects  cause  the  Fermi  level  in  the  “equilibrium”  state 
to  shake  with  time  under  an  ac  electric  field.  As  a  result,  a 
charge-density  fluctuation  in  the  quantum  well  will  modify 
the  Hartree  potential  in  the  surrounding  barrier  region,  and 
thus  greatly  affect  the  sequential  tunneling  of  electrons 
through  the  barrier.  Simultaneously,  the  charge-density  fluc¬ 
tuation  also  modifies  the  electronic  states  in  the  quantum 
well  within  the  self-consistent  Hartree  model.15 

For  the  quantum-well  sample  considered  in  this  study,  the 
second-subband  edge  is  83.4  meV  above  the  first-subband 
edge  (see  Table  I),  while  the  Fermi  level  is  only  14.2  meV 
above  the  first-subband  edge  (see  Table  II).  As  a  result,  the 
second  subband  is  completely  unpopulated  at  temperatures 
below  40  K.  Therefore,  we  have  neglected  the  tunneling 
contribution  from  the  unpopulated  second  subband.  The 
well-known  negative  differential  conductance  (NDC)  phe¬ 
nomenon  can  be  seen  if  the  second  subband  in  the  quantum 
well  is  brought  into  consideration  for  electron  tunneling.  For 
the  multisubband  case,  the  NDC  occurs  at  a  field  strength 
where  the  first-subband  edge  in  the  preceding  well  is  aligned 
with  the  second-subband  edge  in  the  next  tilted  well  (there  is 
a  very  narrow  spectral  density  for  both  quantum  wells  be- 

TABLE  I.  Parameters  of  GaAs/AlxGaj  _x As  MQW  sample  used 
for  numerical  calculations  with  well  depth  Vo,  well  width  LW) 
barrier  thickness  L3 ,  electron  density  n2 d,  well  (barrier)  relative 
dielectric  constant  ew  (eB),  and  well  (barrier)  electron  effective 
mass  mw  (mB)  with  mt  being  the  free-electron  mass. 
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TABLE  n.  Calculated  parameters  of  GaAs/AlxGa!^AsMQW 

sample  used  for  numerical  calculations,  including  average  electron 
effective  mass  m*  with  mt  being  the  tree-electron  mass,  zero-field 
ground-subband  edge  E ,  second-subband  edge  ,  and  chemi¬ 
cal  potential  measured  from  . 

m*  _ (meV) _ (meV)  •  jxQ  (meV) 

0.07  44.5  127.9  14.2 
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lishment  of  the  connection  between  the  current  quantum 
theory  and  the  previous  current- surge  model.  Numerical  re¬ 
sults  and  discussions  are  given  in  Sec.  VI  for  the  changes  of 
current,  drift  velocity,  and  density  as  a  function  of  time.  The 
electron  distribution  functions  at  different  time's  for  various 
temperatures,  amplitudes  of  ac  electric  field,  and  doping  pro¬ 
files  are  also  shown  and  compared.  The  paper  is  concluded  in 
Sec.  VII  with  some  remarks. 


cause  of  the  very  thick  barrier  between  them).  On  the  other 
hand,  the  NDC  phenomenon  also  occurs  in  quantum  wells 
with  a  single  subband.  This  is  due  to  a  lesser  -overlap  be¬ 
tween  the  quantum-well  quasiparticle  spectral  functions 
(whose  width  depends  not  only  on  the  interwell  coupling  but 
also  on  the  disorder  self-energy)  as  the  applied  electric  field 
increases  when  the  Fermi  energy  is  not  too  far  from  the  top 
of  the  barrier.  For  the  sample  with  barrier  thickness  LB 
=  300  A  considered  in  this  paper,  the  required  field  strength 
for  the  multisubband  NDC  phenomenon  is  27:8  kV/cm. 
However,  the  maximum  field  strength  employed  in  this  study 
is  only  1  kV/cm.  Consequently,  we  have  only  included  the 
first  subband  and  neglected  the  NDC  effect  in  our  model 
where  the  Fermi  energy  is  well  below  the  top  of  the  barrier. 
Moreover,  the  field-domain  effect  in  an  MQW  system  is  ex¬ 
pected  to  be  very  small  under  low  electric  fields  for  coherent- 
tunneling  cases  or  below  40  K  for  sequential-tunneling  cases 
and  is  neglected  in  this  paper  since  it  becomes  significant 
only  for  a  large-tunneling  current.  The  thick-barrier-layer 
sample  used  in  this  study  is  to  limit  the  dark  sequential- 
tunneling  current  to  an  extremely  low  amount,  which  ensures 
a  very  high  detectivity  for  quantum-well  infrared  photodetec¬ 
tors  operating  at  a  low-temperature  and/or  a  low-photon 
background."16  When  the  sequential  tunneling  is  low,  the  im- 
punty  or  defect  channels  within  the  barrier  would  play  a 
role.17  However,  this  only  modifies  the  resistance  of  the 
sample  for  sequential  tunneling  of  electrons.  The  nonadia- 
batic  effects  discussed  in  this  paper  for  electron  tunneling 
remain  the  same.  The  usual  self-consistent  Hartree  model  is 
based  on  the  known  equilibrium  (Fenni-Dirac)  distribution 
function  of  electrons,  which  can  be  applied  to  find  electron 
wave  functions  and  energy  levels  simultaneously  in  quantum 
wells.  The  main  result  of  this  paper  is  the  derivation  of  a 
dynamical  differential  equation  for  nonadiabatic  electrons 
under  an  ac  electric  field  in  quantum  wells  which  is  then 
used  to  find  the  electron  distribution  function.  This  dynami¬ 
cal  equation  can  be  coupled  to  the  self-consistent  Hartree 
model  to  solve  for  electron  wave  functions,  energy  levels, 
and  nonadiabatic  distributions  at  the  same  time. 

The  paper  is  organized  as  follows.  In  Sec.  E,  we  introduce 
a  shifted  Fermi-Dirac  model18  for  local  fluctuations  of  elec¬ 
tron  kinetic  energy  and  charge  density  in  the  quantum  well. 
Section  IE  is  used  to  establish  a  unified  theory  for  both  co¬ 
herent  and  sequential  tunneling  of  electrons  in  quantum-well 
and  superlattice  systems.  The  previous  current-surge  model 
is  briefly  reviewed  in  Sec.  IV.  Section  V  is  devoted  to  the 
derivation  of  a  differential  equation  based  on  the  self- 
consistent  Hartree  model  with  the  inclusion  of  nonadiabatic 
effects  on  the  electron  distribution  function,  and  to  the  estab- 


n.  SHIFTED  FERMI-DIRAC  MODEL 

By  using  a  shifted  Fermi-Dirac  model,18  the  nonequilib¬ 
rium  electron  distribution  function  can  be  written  as 


/^(k)=4»(£lk+Ak|+A^),  (1) 


where  f^°(Ek)  is  the  Fermi  function  and  is  given  by 


/?(£*)= 


1  +exp 


Ato(7I2D5n 


k*T 


(2) 


Here,  E  k  is  the  electron  kinetic  energy  and  T  is  the  tempera¬ 
ture.  /io(n2D.n  is  the  chemical  potential  of  the  equilibrium 
electron  gas  relative  to  the  edge  of  the  ground  subband  in  the 
quantum  well  and  is  determined  by  the  electron  density  n2D 
at  T.  In  Eq.  (1),  A Ek  represents  the  local  fluctuation  of  elec¬ 
tron  kinetic  energy  for  electron  state  |k).  Using  the  accelera¬ 
tion  theorem  for  the  momentum  drift  Ak  introduced  in  Eq. 
(1)  under  an  applied  electric  field  £b(t),  we  obtain  the  fol¬ 
lowing  generalized  Boltzmann’s  equation  associated  with  the 
shifted  Fermi-Dirac  model  in  Eq.  (1): 


#*>(k) 

7t 


<?P0(k) 

dEk 


+ 


dLEk  df* o(k) 
dt  dEk 


dt 


coll  > 


(3) 


where  E^t)  £dc^~  £acsin(f2/)  is  the  time-dependent  electric 
field  with  frequency  El  =  2  7rITp9  time  period  Tp ,  dc  ampli¬ 
tude  £dc ,  and  ac  amplitude  E^ .  The  term  on  the  right-hand 
side  of  Eq.  (3)  represents  collision  contributions.  In  the  limit 
of  Dre<a  with  re  being  the  energy-relaxation  time  of  elec¬ 
trons  in  the  quantum  well,  we  have  (d&Ek/dt)t=0;  but 
( dAk/dt)t=eEh(t)!fi ,  where  (♦••),  defines  a  time  average  • 
over  multiple  periods  of  Tp .  In  this  case,  only  the 
momentum-drift  phenomenon  occurs.  On  the  other  hand,  we 
find  (dkkJdt)t=  0  in  the  limit. of  Cirp>  1 ,  with  t  being  the 
momentum-relaxation  time  of  electrons  in  the  quantum  well. 
Under  this  condition,  only  the  energy-drift  phenomenon  oc¬ 
curs  with  {d&Ek/dt)t=  cr(n)£^/4,  where  <r(fl)  is  the  ac 
conductivity  of  electrons  in  the  quantum  well.19 

For  the  shifted  Fermi-Dirac  model  in  Eq.  (1),  there  exists 
a  local  charge-density  fluctuation  for  each  electron  state  Ik) 
defined  by  9 
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yLr»(k) -/?<£*)] 


2e 


ldAEk\ 

\  ^  1 J 

[  dEk  \ 

(4) 


where  V  is  the  volume  of  the  sample.  The  local  charge- 
density  fluctuation  is  a  result  of  the  change  of  the  electron 
distribution  in  energy  space  with  respect  to  the  equilibrium 
State  even  when  n2 d  is  a  constant. 

HI.  COHERENT-  AND  SEQUENTIAL-TUNNELING 
MODELS 

For  an  MQW  system  with  thick-barrier  layers,  the  adia¬ 
batic  sequential-tunneling  current  density  flowing  in  the  z 
direction  (growth  direction  and  perpendicular  to  the 
quantum-well  layers)  is  found  to  be20 


(5) 


where  LB  is  the  thickness  of  the  barrier  between  two  adjacent 
quantum  wells,  v\  is  the  group  velocity  of  quasibound-  or 
continuum-state  electrons  in  the  z  direction,  and  T[Ek,8h]  is  ■ 
the  quantum-mechanical  transmission  of  electrons  through 
the  biased  barrier.  If  Qrt>l,  with  rt  being  the  electron 
sequential-tunneling  time,  T{Eki8{\  has  to  be  found  by 
solving  a  time-dependent  Schrodinger  equation.  Otherwise, 
T[Ek,8h]  can  be  calculated  from  a  static  Schrodinger  equa¬ 
tion  at  each  time  t  if  flrt<^l.  We  will  be  only  interested  in 
the  latter  case  with  1  hereafter. 

In  the  limit  of  small  barrier  thickness  and  weak  field,  i.e., 
e£bLB<Eki  Eq.  (5)  yields  a  coherent- tunneling  current  that 
takes  the  same  form  as  that  obtained  using  the  regular  Bolt¬ 
zmann’s  equation  (dAEk/dt=0)  under  the  relaxation-time 
approximation: 


42  (vl)2rf 


«?/*>(£*) 


BE, 


(6) 


where  we  have  set  T[_Ek,8^vzkrp!LB^  1  for  scattering- 
limited  miniband-state  electrons  at  very  low  electric  field 
(with  mean  free  path  v\rp  smaller  than  LB),  and  rp  is  the 
momentum-relaxation  time  of  electrons.  Consequently,  the 
conductance  That  is  proportional  to  /Ao(t)/£b  becomes  inde¬ 
pendent  of  £b  in  this  situation. 

If  we  replace  f^°(Ek)  to  the  leading-order  approximation 
by  the  equilibrium  value  f%°(Ek)  for  faster  electron  energy 
relaxation  processes  due  to  inelastic  scattering  of  electrons  as 
compared  to  the  electron  sequential  tunneling,  and  replace 
the  electron  group  velocity  v\  by  a  drift  velocity  ud[£b]  (a 
statistically  averaged  group  velocity)  of  electrons  in  a  bulk 
material,  Eq.  (5)  reduces  to  Levine’s  sequential-tunneling 
model20 


^[4]=^ yvlSb]2  T[EkAJ^\Ek) 

~r0\Ek+e  4Lb)], 


(7) 


where  v^8d  =  (eTvlm*)8h,  the  momentum-relaxation  time 
Tp  is  given  by 


(8) 


m*  is  the  effective  mass  of  electrons,  vs  is  the  electron  satu¬ 
ration  velocity,  and  £s  is  the  saturation  electric  field.  In  Eq. 
(7),  ev £5 b]-  can  be  equivalently  viewed  as  a  three- 

dimensional  tunneling-electron  density  that  depends  on  £b, 
r,  and  n2D .  Obviously,  the  conductance  that  is  proportional 
to  Jfl0[8h]/Sh  becomes  dependent  on  8h  in  this  situation. 

IV.  CURRENT-  SURGE  MODEL 

From  now  on,  we  limit  ourselves  to  an  electrical-quantum 
limit  where  only  the  ground  subband  of  the  narrow  quantum 
well  is  occupied  by  electrons  at  low  temperatures  and  low 
electron  densities.  The  electron  kinetic  energy  of  the  ground 
subband  (measured  from  the  edge)  is  given  by  Ek 
In  the  current-surge  model,10’1344  we  assume 
that  A Ek  is  associated  with  the  fluctuation  of  the  chemical 
potential  of  electrons  in  the  quantum  well  (independent  of 
individual  electron  state),  instead  of  the  local  fluctuation  of 
electron  kinetic  energy  for  each  electron  state.  By  writing 
“  A<u=^,0(w2d  At(0  for  the  global  chemical- 
potential  fluctuation,  where  Jt)  and  fiio(n2I>,T)  are,  respec¬ 
tively,  the  transient  chemical  potential  for  electron  density 
ne(0  that  for  an  equilibrium  electron  gas  in  quantum 
wells,  we  get 


dAEk ^  dA Ek  ^  djx  dnt 
dt  dt  ■  dll,,  dt 


(9) 


We  further  introduce  a  spatially  averaged  space-charge  field 
8Jt)  which  is  defined  by10,1344 


A Ek  e8Jt)LB, 


(10) 


where  8  Jit)  measures  the  reduction  of  the  electron  chemical 
potential  in  quantum  wells.  If  we  use  Levine’s  sequential- 
tunneling  model  in  Eq.  (7),  we  find  the  change  in  the  current 
density  due  to  the  existence  of  this,  space-charge  field  Sjt), 


(11) 


where  JMo[£b]  has  been  given  in  Eq.  (7).  In  Eq.  (11),  the  first 
term  can  be  viewed  as  an  equivalent  capture  current  flowing 
into  the  quantum  well,  while  the  second  term  can  be  re¬ 
garded  as  a  sequential-tunneling  current  flowing  out  of  the 
quantum  well. 

For  a  quantum  well,  the  electron  density  will  be  constant 
if  the  conduction  currents  flowing  in  and  out  of  the  well  are 
equal.  The  variation  of  the  charge  density  in  the  well  is  cre- 
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ated  by  an  imbalance  in  conduction  currents.  The  charge- 
current  conservation  law  for  density  fluctuation  8p(t)  re¬ 
quires 

Spk{t)=AKJM,  (12) 

where  A  is  the  cross-sectional  area  of  the  MQW  sample.  The 
left-hand  side  of  Eq.  (12)  represents  the  charge  increase  in¬ 
side  the  well,  while  the  right-hand  side  of  the  equation  stands 
for  the  net  increase  in  charges  due  to  a  nonadiabatic  change 
in  the  current  flowing  into  the  quantum  Well.  Combining 
Eqs.  (4),  (9),  (10),  and  (12),  we  finally  arrive  at  the  following 
equation  derived  previously  as  the  current-surge  model:10 
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doping  for  the  single  quantum  well,  n&(z,t) 
=  \4>\(z,t)\2nt(t)  is  the  density  function,  and 

Sp(t)Lvj  f+=e 

«e(0  =  «2D+ - =  ^2D+P2D  dE8f(E,t), 

e  Jo 

(16) 

where  p2D=(mw/'77*^2)  is  the  density  of  states  for  two- 
dimensional  electrons  in  the  quantum  well  and  8f(E,t)  rep¬ 
resents  the  local  fluctuation  of  the  electron  distribution  func¬ 
tion  in  energy  space.  Here,  the  number  of  electrons  in  the 
quantum  well  is  not  a  constant  due  to  the  nonadiabatic  cur¬ 
rent  flowing.  Moreover,  we  find  from  Eqs.  (4),  (9),  and  (16) 
that 


d  d 

L*CW'dt^t}~L*C  Qw^b(0  ”  AJna(0>  (13) 

where  the  quantum-well  capacitance  per  unit  area  is  CQW 
=  (m*e2/TTh2)f%°(0).  Here,  we  have  employed  in  Eq.  (13) 
the  fact  that  (dp,/ dn^)(dn&/dt)  =  eLB[dSh(t)/dt ]  for  the  ca¬ 
pacitance  coupling  of  the  quantum  well  to  an  external  ac 
electric  field.10,13,14  The  fast  inelastic  scattering  in  quantum 
wells  ensures  that  electrons  are  in  an  “equilibrium”  state. 
However,  Eq.  (13)  causes  a  shaking  Fermi  level  for  the  equi¬ 
librium  state  on  a  macroscopic  time  scale. 


^wdt  ~~  ef>2D  J0  ~st 

o  d  f +cc 

+  *%P2D^4(0  Jo  dE 

(17) 

Applying  Eq.  (12)  and  using  Eqs.  (11)  and  (17),  we  find  the 
following  integral  equation  for  8f(Eyt)  by  using  Levine’s 
model  in  Eq.  (7) 


<°(£) 

dE 


V.  NONADIABATIC  EFFECTS  IN  A  SELF-CONSISTENT 
HARTREE  MODEL 

As  mentioned  in  the  Introduction,  electrons  in  quantum 
wells  see  only  the  instantaneous  ac  electric  field  during  their 
sequential-tunneling  process  if  L  In  this  case,  the 
ground-state  electron  wave  function  4>i(z,t)  inside  the  quan¬ 
tum  well  within  the  self-consistent  Hartree  model  is  deter¬ 
mined  by21 


h2  d  I  1  d  \ 


+vH (*,<) 


(14) 


where  E}(t)  is  the  time- dependent  ground-subband  edge,  the 
electron  effective  mass  m*(z)  takes  mw  in  the  well  and  mB 
in  the  barrier,  and  Z7qW(z)  is  zero  inside  the  well  but  F0 
outside  the  well.  For  the  adiabatic  state,  we  have  p(t) 
otherwise  8p(t)i=0  for  the  cases  with  nona¬ 
diabatic  effects.  The  Hartree  potential  VH(z,t)  .in  Eq.  (14) 
can  be  found  from  the  Poisson  equation 


d  , 

eP2Djo  dE—Sf(E,l)  +  e-LBp2I) 


dm 

dt 


"  <°(E)' 

/  eP2T)\ 

dE 

\  / 

W41+ 


xf 

X  ^dET{E+Ex  ,4;Vh]L f0\E)  +  8f{E,i) 
-fQ\E + e£hLB) -  Sf{E + eS ^ ,f)] 

+  ( v£Sd \+J'dET[E+Ef>  A'M^^iE) 


-/*>(£+ e4LB)]  =  0. 


(18) 


In  Eq.  (18),  VH(z,t)  and  E^t)  are  written  simply  as  VH  and 
Ex.  The  adiabatic  quantities  VBJ(z,t)  and  E[0>(t)  can  be 
obtained  by  simply  setting  Sp(t)  =  0  in  Eq.  (15)  and 
VH  (z,t)  =  V%\z  ,t)  in  Eq.  (14).  Moreover,  the  fluctuation  of 
the  drift  velocity  8vd[8f\  introduced  in  Eq.  (18)  is  calcu¬ 
lated  to  be 


e 2 

=  7'[^D(z)-«eC  Z,t)]>  (15) 

e0 

where  donors  are  assumed  completely  ionized,  and  the  rela¬ 
tive  dielectric  constant  er(z)  takes  ew  in  the  well  and  eB  in 
the  barrier.  NB(z)  in  Eq.  (15)  is  the  static  profile  of  donor 


/  P2d\  f4*03  f2E 

8vl8f]  =  -\—\  dE8f(Eyt )V— •  (19) 

\n2D/  JO  y 

Finally,  Eq.  (18)  leads  us  to  the  dynamical  differential  equa¬ 
tion  for  8f(E,t), 


d 

dz 


er(z)  ^VH(z,0 
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jtdf(E,t)-eLh 


dm  <°(-g) 

dt  BE 


Aw 


Cqw—  e2  pm 


<“(£)' 

dE  ’ 


-X{vim8vl8f]}T[E+E,  ,4;Vh]  [/*>(£) 

+  Sf(E,t)  -fZ°(E  +  'em)  ~  SfiE+eS^  ,t)] 

+  ^vl^nE+E?\£h;V^[fZ°(E) 

-f0\E+e  4Lb)]=0,  (20) 

where  the  initial  condition  is  chosen  to  be  8f{E,  t)~  0  at  t 
=  0  if  the  ac  electric  field  is  applied  to  the  system  after  t 
=  0.  8f{E,t)  has  a  lower  bound  that  is  set  by  the  condition 

mo  +4°(^)=o. 

For  small  A/z,,  the  first  term  in  Eq.  (20)  can  be  approxi¬ 
mated  to  the  leading  order  by 


^[4]=(^)vd[4]/o+W[E+^)>4;v^] 

X[C(£)-/T(£+e4LB)], 

T[E + £<°> ,  4 ;  Vg>  -eEmtLBi^T[E+ E™  ,  4+  4g  V^] , 
and  ud[£b]  +  8v  ^8f]^v  8^  if  we  set  rp 

^LBf(v¥f 2)  with  uF  being  the  electron  group  velocity  at  the 
Fermi  level.  The  space-charge  field  8m(t)  introduced  by 
A/U,=  “  e^na(0-^B  can  be  calculated  from 


at 


8f(E,t)~ 


d&jj, 

dt 


df0\E) 

BE 


(21) 


Similarly,  a  part  of  the  third  term  in  Eq.  (20)  can  be  approxi¬ 
mated  as 


which  becomes  positive  if  8p(i)< 0. 

The  quantum-mechanical  transmission  coefficient  T[E 
+E}  ,4;Vh]  used  to  evaluate  the  tunneling  current  in  Eq. 
(7)  can  be  found  from  the  following  backward  iteration8  at 
each  time  t 


TiE+E,  AlVui[f^\E)  +  8f(E,t) 

-/£°(£+e4LB)-  <?/(£+ e4LB,f)] 


*/-!«=  2+F[C/B-e4(0(/-l)A  +  Pf(f)-£' 


T[E: + > ,  4 ;  V^0) + ( <5V^0)/  <5n2D)p2DA/i] 

X[^“+A4E)-^+A4£+e5b£B)];  (22) 


where  (<?V^0)/5«2d)  =  (e2/2e0ew?TF)  and  qrF 

=  (e2/2e0ew)p2D  in  the  Thomas-Fermi  model.21  Byrecalling 
A^=-e£na(Y)Z,Bj  Eq.  (20)  results  in  the  current- surge 
model  in  Eq.  (13),  where 


-e  m 


(24) 


for where  ^(f)=  ^(0=^,0, 

Ea—h  /2tmbA  ,•  A=Z.B/7VB  and  jVb  is  the  number  of  slabs 
(thickness  A)  within  the  barrier  layer.  Here,  I/®=0  for  j 
=  0  andy'=AB+l.  Otherwise,  Uj=V0.  The  ending  bound- 
ary  condition  of  Eq.  (24)  produces 


lM 

^ArB(0 


exp(zi\fB£') 

■C1  —  C'S-b+i  +  «4CO-2VbA  — F^b+1(0D>^bh-iW 


(25) 


where  k'  —  (A/ft)  y/2ms{E + e£b(t)L^) .  From  the  solution  of 
Eq.  (24)  we  find  the  quantum-mechanical  transmission  of 
electrons  from 


a 

A 

~(i/2  k){4>2{t)~4>o{t)). 

(27) 


T[E+EU  4;Vh]  = 


_1  _  jE+e£b(t)LB 
\S\2\  E 


(26) 


where  |5,|2  =  [|a|2  +  |fi|2+2Re(a6*)]/4.  The  solution  of  Eq. 
(24)  ensures  that  the  transmission  coefficient  in  Eq.  (26)  de¬ 
pends  on  the  barrier  thickness  and  height  in  an  exponential 
way.  Here,  the  two  complex  numbers  a  and  b  are  defined  by 
the  starting  boundary  condition  of  Eq.  (24) 


with  k=  (A IK)  Jlm^E. 

VT.  NUMERICAL  RESULTS  AND  DISCUSSION 

We  choose  a  GaAs/ AlxGa^  _xAs  MQW  sample  for  nu¬ 
merical  calculations.  Some  sample  parameters  can  be  found 
in  Tables  I  and  n.  Other  parameters  include  4c 
=  0.05kV/cm,  us=2X106  cm/sec,  4=2  kV/cm,  and  Tp 
—  4  sec.  Different  doping  profiles  have  been  considered,  in- 
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FIG.  1.  Calculated  time  dependence  of  electron  sequential- 
tunneling  currents  I(t)  =  J[JfLa[£h]  +  (left  scale)  for  nona- 

diabatic  [AV^O^O,  solid  curve]  and  adiabatic  [AJm(t)=0,  '■ 
dashed  curve],  and  an  applied  ac  electric  field  £b(f)  (right  scale)  for 
T—  40  K  and  £^—1  kV/cm  with  uniform  doping  inside  GaAs 
quantum  wells. 

eluding  uniform  well  doping,  center  8  doping,  and  edge  8 
doping.  The  temperature  T  and  the  amplitude  of  the  ac  elec¬ 
tric  field,  £ac,  will  be  given  in  the  figure  captions. 

Figure  1  presents  the  calculated  nonadiabatic  electron 
sequential-tunneling  current  /na(0  =  [^°[4]  +  AJna(0]-A 
(solid  curve)  and  the  adiabatic  sequential-tunneling  current 
(dashed  curve)  as  a  function  of  time  t  (left 
scale)  for  the  MQW  sample  at  T=40  K,  1  kV/cm,  and 
with  uniform  doping  inside  the  GaAs  quantum  wells.  For 
comparison,  the  applied  ac  electric  field  £b(i)  (right  scale)  is 
also  plotted  in  the  same  figure.  When  £b(/)  approaches  its 
maximum  (i.e.,  t=\  sec),  we  find  a  small  enhancement  in 
with  respect  to  7a(f)  and  the  saturation  of  IrJt)  due  to 
the  large  reduction  in  electron  density  inside  the  quantum 
wells.  On  the  other  hand,  we  find  a  large  enhancement  in 
/na(0  due  to  the  large  increase  of  electron  density  inside  the 
quantum  wells  when  £b(t)  approaches  its  minimum  (i.e.,  t 


=3  sec).  These  features  are  a  result  of  the  induced  space- 
charge  field  . 

The  reduction  of  electron  -population  around  E 
=  yU0(w2D>^)  can  be  described  by  the  space-charge  field 
£na(r)  defined  in  Eq.  (23).  We  display  EjJj)  (solid  curve,  left 
scale)  in  Fig.  2(a),  along  with  £b(f)  (dashed  curve,  right 
scale)  as  functions  of  t  for  uniform  doping.  From  the  figure 
we  see  that  £na(r)  and  £b(t)  are  nearly  in  phase  with  each 
other,  except  for  a  slight  phase  shift.  This  is  a  direct  result  of 
oscillations  in  the  change  of  the  charge  density  8p(t)  in  the 
quantum  well,  as  shown  in  Fig.  2(b),  where  both  8p(t)  (solid 
curve,  left  scale)  and  8v£8f]  (dashed  curve,  right  scale)  are 
plotted  as  functions  of  t.  Since  SjJJ)  describes  the  reduction 
of  charge  density  in  the  quantum  well,  we  expect  8p(t)  to  be 
nearly  out  of  phase  with  .£^(0  or  £b(t),  as  can  be  seen  from 
Figs.  2(a)  and  2(b).  The  situations  with  Sp(t)< 0  and 
8p(t)>0  indicate  electrons  moving  out  of  and  into  the  quan¬ 
tum  well,  respectively.  Moreover,  8v£8f]  will  be  in  phase 
with  £^{t)  since  it  is  proportional  to  —  8f{E,t)  that  itself  is 
proportional  to  £na(  t). 

Figure  3  displays  the  calculated  adiabatic  Hartree  poten¬ 
tial  [in  panel  (a)]  and  the  change  of  Hartree  potential  [in 
panel  (b)]  in  the  nonadiabatic  state  from  the  Poisson  equation 
(15)  as  functions  of  position  z  for  different  doping  profiles  at 
t!Tp  =  0.25.  Here,  T—  40  K,  £^=1  kV/cm,  and  the  quantum 
well  sits  in  the  range  of  300  A=£z^380  A.  From  Fig.  3(a) 
we  find  that  the  absolute  value  of  the  adiabatic  Hartree  po¬ 
tential  becomes  smallest  for  the  uniform-doping  case.  The 
center  8  doping  in  the  quantum  well  causes  the  conduction 
band  edge  to  bend  down  at  the  well  center,  while  the  edge  8 
doping  makes  the  conduction  band  edge  bend  up  there,  as 
shown  in  Fig.  3(a).  With  the  total  potential  seen  by  the  elec¬ 
trons  being  the  sum  of  the  adiabatic  Hartree  potential 
F^O,0  plus  the  change  A plus  the  quantum-well 
potential  f/QW(z),  the  out-of-phase  nature  of  Figs.  3(a)  and 
3(b)  will  result  in  the  band  bending  seen  in  Fig.  3(a)  being 
substantially  suppressed  by  the  nonadiabatic  effects  in  Fig. 
3(b).  However,  the  nonadiabatic  effects  with  edge  8  doping 
produces  two  positive  spikes  [solid  curve  in  Fig.  3(b)]  at  the 


t  ( sec ) 

FIG.  2.  Time  dependence  of  a  calculated  space-charge  field  £naW  (solid  curve,  left  scale)  and  an' applied  ac  electric  field  £b(t)  (dashed 
curve,  right  scale)  in  panel  (a),  as  well  as  time  dependence  of  calculated  charge-density  fluctuation  8p(t)  (solid  curve,  left  scale)  and  change 
of  drift  velocity  (5U{j[  <5/1  (dashed  curve,  right  scale)  in  panel  (b)  for  T—  40  K  and  £ac=  1  kV/cm  with  uniform  doping  inside  quantum  wells. 
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FIG.  3.  Calculated  position  dependence  of  adiabatic  Hartree  potentials  V^\z,ty  in ‘panel  (a)  and  change  of  nonadiabatic  Hartree 
potentials  A  in  panel  (b)  at  t/Tp= 0.25,  with  F=40  K  and  1  kV/cm  for  edge  S  doping  (solid  carves),  center  £  doping  (dashed 

curves),  and  uniform  doping  (dotted  curves)  in  the  quantum  well. 


edges  of  the  quantum  well  and  thus  will  reduce  the  electron 
sequential-tunneling  current. 

Figure  4  indicates  the  effects  of  an  ac  electric-field 
strength  £ac  [in  panel  (a)]  and  temperature  T  [in  panel  (b)]  on 
the:  charge-density  fluctuations  Sp(t)  as  a  function  of  t  in  a 
uniformly  doped  quantum  well.  In  Fig.  4(a)  we  find  that 
fluctuations  Sp(t)  increase  with  at  T-  40  K,  with  the 
negative  peak  (electrons  removed,  from  the  quantum  well) 
being  saturated  at  8^=5  kV/cm.  In  Fig.  4(b),  as  T7  increases 
5p{t)  is  enhanced  when  it  is  negative  (electrons  removed 
from  the  quantum  well),  but  reduced  when  it  is  positive 
(electrons  added  to  the  well)  at  8^=  1  kV/cm. 

Figure  4  only  shows  us  the  global  fluctuation  of  the 
charge  density  in  the  quantum  well.  In  order  to  gain  further 
insight  into  the  local  change  in  the  electron  distribution  func¬ 
tion,  we  display  8f{E,t )  in  Fig.  5  at  ?/7^  =  0.25  with  uni¬ 
form  doping  for  different  values  of  £ac  [in  panel  (a)]  and  T . 


t  (sec) 


[in  panel  (b)].  From  Fig.  5(b)  it  is  clear  that  Sf(E,i)  always 
shows  a  negative  minimum  at  Po(n2D,T),  since  it  is  propor¬ 
tional  to  df%°(E)/dE.  Because  the  Fermi  surface  broadens 
with  increasing  T9  we  find  from  Fig.  5(b)  that  the  negative 
minimum  is  partially  suppressed  and  broadened  (solid  curve) 
when  T—  40  K  as  compared  to  that  (dashed  curve)  at  T 
~  20  K.  From  Fig.  5(a)  we  find  that  the  negative  minimum  is 
enhanced  when  8^  is  increased.  The  cusp  (dashed  curve)  in 
Fig.  5(a)  is  a  result  of  zero  occupation  of  electrons  in  a 
specific  state  with  kinetic  energy  E  in  the  ground  subband. 

Figures  6(a)  and  6(b)  present  nonadiabatic  effects  on  the 
Hartree  potentials  in  the  uniformly  doped  quantum  well  at 
r=40  K  and  £ac=  1  kV/cm.  From  Fig.  6(a)  we  find  that  the 
positive  peak  in  the  adiabatic  Hartree  potential  V$\z,t) 
(dashed  curve)  at  the  center  of  the  quantum  well  is  greatly 
suppressed  by  the  nonadiabatic  effects  (solid  curve)  at  tIT 
=  0.25,  leaving  two  positive  spikes  at  the  edges  of  the  quan- 


0  1  2  3  4 


t  ( sec ) 


FIG,  4.  Time  dependence  of  calculated  charge-density  fluctuations  Sp(t)  in  the  uniformly  doped  quantum  well  .  In  panel  (a),  we  set 
T=40  K  with  E^=  1  kV/cm  (solid  curve)  and  £^=  5  kV/cm  (dashed  curve),  in  panel  (b),  we  set  £ac=  1  kV/cm  with  T=  40  K  (solid  curve) 
and  20  K  (dashed  curve). 
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FIG.  5.  Calculated  change  of  nonadiabatic  distribution  functions  i 
weU.  In  panel  (a),  we  set  T=  40  K  with  £K=1  kV/cm  K  (solid  c 
=  1'  kV/cm  with  T=  40  K  (solid  curve)  and  7=  20  K  (dashed  curve). 


0  10  20  30 
E  ( meV ) 

',/)  at  t/Tp=0.25  for  electrons  in  the  uniformly  doped  quantum 
)  and  £ac=  5  kV/cm  (dashed  curve).  In  panel  (b),  we  set  8^ 


turn  well.  Figure  6(b)  shows  the  comparison  between  Hartree 
potentials  when  the  electrons  in  the  quantum  well  are  either 
removed  (t/Tp  =  0.25,  solid  curve)  or  added  (f/r  =0.75, 
dashed  curve).  We  find  from  Fig.  6(b)  that  the  two  positive 
spikes  at  the  edges  of  the  quantum  well  are  suppressed,  but 
two  negative  spikes  are  generated  when  electrons  are  added 
to  the  well. 

Finally,  we  display  in  Fig.  7  f(E,t)  at  ftTp=0.25  (dotted 
curve)  and  0.75  (dashed  curve),  as  well  as  the  equilibrium 
distribution  f%\E)  (solid  curve)  in  panel  (a)  and 
l°gioUna(Ol  as  a  function  of  £b(t)  in  panel  (b).  From  Fig. 
7(a)  we  see  f(E,t)  resembles  the  equilibrium  distribution 
with  a  shaking  Fermi  level  with  time  (shaking  up  at 
t/Tp— 0.75  and  shaking  down  at  t/Tp=0.2S).  Compared 
with  the  adiabatic  electron  sequential-tunneling  current  [thin 
solid  curve  with  AJna(f)  =  0]  in  Fig.  7(b),  the  symmetry  of 
i°gioP*na(Ol  with  respect  to  the  positive  (electrons  being  re¬ 


moved)  and  negative  (electrons  being  added)  extreme  values 
of  £b(t)  is  broken  in  the  case  with  nonadiabatic  effects  (thick 
solid  curve).  A  small  offset10-13-14  of  log10|/na(t)|  with  respect 
to  £b(0  =  0  can  be  seen  by  comparing  thick  and  thin  solid 
curves. 


vn.  CONCLUSIONS  AND  REMARKS 

In  conclusion,  we  have  derived  a  dynamical  differential 
equation  for  the  nonadiabatic  electron  distribution  function 
with  sequential-tunneling  current  flowing  through  an  MQW 
system.  Using  this  equation,  we  generalized  the  self- 
consistent  Hartree  model  for  the  calculation  of  the  electronic 
states  with  the  inclusion  of  nonadiabatic  effects  in  a  quantum 
well.  We  have  also  studied  the  effects  of  different  doping 
profiles,  temperatures,  and  amplitudes  of  applied  ac  electric 
field  on  the  nonadiabatic  electron  sequential  tunneling  Fi- 


z  (A) 


^FIG.  6.  Calculated  position  dependence  of  Hartree  potentials  for  nonadiabatic  (solid  curve)  and  adiabatic  (dashed  curve)  in  panel  (a)  at- 
7-40  K  and  £*.=  1  kV/cm  with  uniform  doping  inside  the  quantum  well,  and  nonadiabatic  Hartree  potentials  FH(z, t)  in  panel  (b)  for 
t/Tp  —  Q.25  (solid  curve)  and.  t! Tp  =  0.75  (dashed  curve). 
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E  (meV)  eb(t)  (kV/cm) 


^  nonadiabatic  electron  distribution  functions  f{E,t)  in  panel  (a)  and  logarithm  of  absolute  value  of  nonadiabatic 

C  q= ,  ,  v/  m®kng  C™ef  “[Sbl  +  A/^O]  as  a  function  of  an  applied  ac  electric  field  4(r)  in  panel  (b)  for  T=40  K  and 

“  ..  4  ’  W,e  p  °t/(jE,f)  at  tlTP=0:25  (dotted  curve)  and  tlTp= 0,75  (dashed  curve).  The  equilibrium  electron  distribution 

nTw5  2- 1  T?}  18  ShOWn  f0r  *e  COmParison-  1*  P^61  Cb),  the  currents  calculated  from  adiabatic  (thin  curve)  and 

nonadiabatic  (thick  curve)  electron  sequential  tunneling  are  compared  with  each  other. 


nally,  we  have  connected  the  present  quantum-statistical 
theory  to  the  previously  proposed  current-surge  model  with  a 
leading-order  approximation. 

In  this  paper,  only  the  self-consistent  Hartree  model  is 
employed.  The  exchange  interaction  between  electrons  and 
the  field-domain  effect  are  expected  to  be  very  small22  at  T 
=40  K  and  £‘ai=  1  kV/cm  and  have  been  neglected. 

The  time  scale  for  observing  the  nonadiabatic  effects  re¬ 
quires 


re<Tt<t<RtC(yffA,  2  t r/Xl, 

where  d&^  1  is  the  differential  tun¬ 

neling  resistance,  depending  on  4  and  T.  Here,  Xlre<8  1  ex¬ 
cludes  the  energy-drift  effect,  leaving  only  the  momentum- 
drift  effect  in  the  system.  Further,  fit, <81  ensures  that  the 
electrons  see  only  an  instantaneous  ac  electric  field  during 
the  sequential-tunneling  process.  Finally,  f<E,CQWyl  en¬ 


sures  the  observation  of  the  nonadiabatic  effects  inside  the 
quantum  well.  Assuming  re=  1  ps  corresponds  to  a  homoge¬ 
neous  energy-level  broadening  of  1  meV,  leading  to  XI 
^10  Hz  from  XI  te<8 1 .  Therefore,  only  a  momentum  drift 
exists  for  low  ac  frequency  Q~ 1  Hz.  The  tunneling  time  rt 
can  be  estimated  from  r(~e/[^[4]^].  For  a  superlattice, 
we  take  J^[£^\A=\  /r A,  leading  to  r,=  0.1  ps  and  XI 
^  10  Hz  from  Hr^l.  For  an  MQW  system,  we  take 
-N1°[4]-4=  10  pA,  leading  to  rt=  10  ns  and  Xl<8108  Hz. 
This  justifies  the  calculation  of  the  quantum-mechanical 
transmission  of  electrons  through  a  biased  barrier  using  the 
time-independent  Schrodinger  equation  in  an  MQW  system 
with  XI  1  Hz.  Difficulties  in  observing  the  nonadiabatic  ef¬ 
fects  may  come  from  the  small  quantum-well  capacitance 
cqw-4~  10  pF  in  the  requirement  t<R,CQWA  For  a  super- 
lattice,  we  take  i?,=  104  ohm,  and  then  f<10~7  sec  is  re¬ 
quired  (impossible  to  observe  with  XI ~1  Hz).  For  an  MQW 
system,  on  the  other  hand,  we  take  R,=  1011  ohm,  which 
implies  t<\  sec  (very  easy  to  observe  with  XI ~1  Hz). 
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